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Discovery of Seven Novel Mammalian and Avian Coronaviruses in the 
Genus Deltacoronavirus Supports Bat Coronaviruses as the Gene 
Source of Alphacoronavirus and Betacoronavirus and Avian 
Coronaviruses as the Gene Source of Gammacoronavirus and 
Deltacoronavirus 
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Recently, we reported the discovery of three novel coronaviruses, bulbul coronavirus HKU11, thrush coronavirus HKU12, and 
munia coronavirus HKU13, which were identified as representatives of a novel genus, Deltacoronavirus, in the subfamily Coro- 
navirinae. In this territory-wide molecular epidemiology study involving 3,137 mammals and 3,298 birds, we discovered seven 
additional novel deltacoronaviruses in pigs and birds, which we named porcine coronavirus HKU15, white-eye coronavirus 
HKUI16, sparrow coronavirus HKU17, magpie robin coronavirus HKU18, night heron coronavirus HKU19, wigeon coronavirus 
HKU20, and common moorhen coronavirus HKU21. Complete genome sequencing and comparative genome analysis showed 
that the avian and mammalian deltacoronaviruses have similar genome characteristics and structures. They all have rela- 
tively small genomes (25.421 to 26.674 kb), the smallest among all coronaviruses. They all have a single papain-like pro- 
tease domain in the nsp3 gene; an accessory gene, NS6 open reading frame (ORF), located between the M and N genes; and 
a variable number of accessory genes (up to four) downstream of the N gene. Moreover, they all have the same putative 
transcription regulatory sequence of ACACCA. Molecular clock analysis showed that the most recent common ancestor of 
all coronaviruses was estimated at approximately 8100 BC, and those of Alphacoronavirus, Betacoronavirus, Gammacoro- 
navirus, and Deltacoronavirus were at approximately 2400 BC, 3300 BC, 2800 BC, and 3000 BC, respectively. From our 
studies, it appears that bats and birds, the warm blooded flying vertebrates, are ideal hosts for the coronavirus gene source, 
bats for Alphacoronavirus and Betacoronavirus and birds for Gammacoronavirus and Deltacoronavirus, to fuel coronavi- 


rus evolution and dissemination. 


oronaviruses (CoVs) are found in a wide variety of animals, in 

which they can cause respiratory, enteric, hepatic, and neuro- 
logical diseases of varying severity. Based on genotypic and sero- 
logical characterization, CoVs were traditionally divided into 
three distinct groups (3, 22, 54). Recently, the Coronavirus Study 
Group of the International Committee for Taxonomy of Viruses 
has proposed three genera, Alphacoronavirus, Betacoronavirus, 
and Gammacoronavirus, to replace the traditional CoV groups 1, 
2, and 3. As a result of the unique mechanism of viral replication, 
CoVs have a high frequency of recombination (22). Their ten- 
dency for recombination and the inherently high mutation rates 
in RNA virus may allow them to adapt to new hosts and ecological 
niches (18, 47). 

The recent severe acute respiratory syndrome (SARS) epidemic, 
the discovery of SARS coronavirus (SARS-CoV), and the identifica- 
tion of SARS-CoV-like viruses from Himalayan palm civets and a 
raccoon dog from wild live markets in China have boosted interest in 
the discovery of novel CoVs in both humans and animals (5, 16, 33, 
36, 39, 40, 46). A novel human CoV (HCoV) of the genus Alpha- 
coronavirus, human coronavirus NL63 (HCoV-NL63), was re- 
ported independently by two groups in 2004 (12, 44). In 2005, we 
also described the discovery, complete genome sequence, clinical 
features, and molecular epidemiology of another novel HCoV, 
human coronavirus HKU1 (HCoV-HKU1), in the genus Beta- 
coronavirus (24, 48, 50). As for animal CoVs, we and others have 
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described the discovery of SARS-CoV-like viruses in horseshoe 
bats in Hong Kong Special Administrative Region (HKSAR) and 
other provinces of China (25, 30). Based on these findings, we 
conducted molecular surveillance studies to examine the diversity 
of CoVs in bats of our locality as well as of the Guangdong prov- 
ince of southern China, where the SARS epidemic originated and 
wet markets and game food restaurants serving bat dishes are 
commonly found. In these studies, at least nine other novel CoVs 
were discovered, including two novel subgroups in Betacoro- 
navirus, subgroups C and D (26, 37, 45, 51). Other groups have 
also conducted molecular surveillance studies in bats and other 
animals, and additional novel CoVs were discovered and com- 
plete genomes sequenced (4, 6, 7, 9, 10, 13-15, 17, 21, 31, 32, 
34, 43, 53). 

Birds are the reservoir of major emerging viruses, most no- 
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TABLE 1 Animals screened and associated CoVs in the present surveillance study 


No. of No. (%) of specimens 

Animal Sample type specimens tested positive for CoV CoV 

Asian leopard cat Rectal swab and tracheal swab 30 0 

Bat Rectal swab 434 0 

Bird* Rectal swab 3,306 35 (1.1%) WECoV HKU16 (n = 3), SpCoV HKU17 (n = 7), 
MRCoV HKUI18 (n = 1), NHCoV HKU19 
(n = 5), WiCoV HKU20 (n = 1), CMCoV 
HKU21 (n = 1), BuCoV HKU11 (n = 10), 
ThCoV HKU12 (n = 1), MunCoV HKU13 
(n = 6) 

Cat Rectal swab and tracheal swab 460 0 

Cattle Rectal swab 47 0 

Chicken Cloacal swab 221 0 

Dog Rectal swab and tracheal swab 462 0 

Human NPA? 1,387 0 

Monkey Rectal swab 235 0 

Pig Rectal swab 169 17 (10.1%) PorCoV HKUI15 

Rodent Rectal swab 389 0 


* No. of birds tested for individual species and their associated CoVs are listed in Table S2 in the supplemental material. 


» NPA, nasopharyngeal aspirate. 


tably, avian influenza viruses (29). Due to their flocking behav- 
ior and abilities to fly over long distances, birds have the po- 
tential to disseminate these emerging viruses efficiently among 
themselves and to other animals and humans. As for CoVs, the 
number of known CoVs in birds is relatively small compared to 
that in bats. Recently, we described the discovery of three novel 
CoVs in three families of birds, named bulbul coronavirus 
HKU11 (BuCoV HKU11), thrush coronavirus HKU12 (ThCoV 
HKU12), and munia coronavirus HKU13 (MunCoV HKU13) 
(49). These three CoVs formed a unique group of CoV, which 
probably represented a novel genus of CoV, Deltacoronavirus (8). 
We hypothesize that there are other previously unrecognized 
CoVs in this novel genus from mammals and other families of 
birds. To test this hypothesis, we carried out a territory-wide 
molecular epidemiology study in 3,137 mammals and 3,519 
birds in HKSAR. Based on the results of comparative genome 
and phylogenetic analysis in the present study, we propose 
seven novel CoVs in Deltacoronavirus. Our model of bats and 
birds as the gene source of the four genera of coronaviruses is 
also discussed. 


MATERIALS AND METHODS 


Animal surveillance and sample collection. All specimens of bats, cats, 
dogs, wild rodents, monkeys, and birds were collected with the assistance 
of the Department of Agriculture, Fisheries and Conservation, HKSAR, 
and those of pigs, cattle, chickens, and street rodents were collected with 
the assistance of the Department of Food, Environmental and Hygiene, 
HKSAR, from various locations in HKSAR over a 53-month period (Feb- 
ruary 2007 to June 2011). All specimens of Asian leopard cats were col- 
lected in the Guangdong province of southern China over an 8-month 
period (August 2010 to March 2011). Tracheal, rectal, and cloacal swabs 
were collected using procedures described previously (47, 49). Nasopha- 
ryngeal aspirates from humans were collected from patients in Queen 
Mary Hospital over a 13-month period (February 2010 to February 2011) 
(24, 47, 50). A total of 7,140 samples from 11 species of bats, 169 pigs, 230 
cats, 231 dogs, 47 cattle, 221 chickens, 389 rodents, 235 monkeys, 1,397 
humans, 15 Asian leopard cats, and 3,298 dead wild birds of 134 different 
species in 38 families had been tested. 

RNA extraction. Viral RNA was extracted from the tracheal, rectal, 
and cloacal swabs and nasopharyngeal aspirates using RNeasy Mini Spin 
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column (Qiagen, Hilden, Germany) (27, 45, 47, 50). The RNA was eluted 
in 50 pl of RNase-free water and was used as the template for reverse 
transcription-PCR (RT-PCR). 

RT-PCR of RdRp gene of CoVs using Deltacoronavirus conserved 
primers and DNA sequencing. Initial CoV screening was performed by 
amplifying a 440-bp fragment of the RNA-dependent RNA polymerase 
(RdRp) gene of CoVs using Deltacoronavirus conserved primers (5’-GTG 
GVTGTMTTAATGCACAGTC-3’ and 5’-TACTGYCTGTTRGTCATRG 
TG-3’) designed by multiple alignments of the nucleotide sequences of 
available RdRp genes of BuCoV HKU11, ThCoV HKU12, and MunCoV 
HKU13 (49). Reverse transcription was performed using the SuperScript 
III kit (Invitrogen, San Diego, CA). The PCR mixture (25 pl) contained 
cDNA, PCR buffer (10 mM Tris-HCl, pH 8.3, 50 mM KCl, 3 mM MgCL, 
and 0.01% gelatin), 200 wM each deoxynucleoside triphosphate (dNTP), 
and 1.0 U Taq polymerase (Applied Biosystems, Foster City, CA). The 
mixtures were amplified with 60 cycles of 94°C for 1 min, 48°C for 1 min, 
and 72°C for 1 min and a final extension at 72°C for 10 min in an auto- 
mated thermal cycler (Applied Biosystems, Foster City, CA). Standard 
precautions were taken to avoid PCR contamination, and no false positive 
was observed in negative controls. 

The PCR products were gel purified using the QIAquick gel extraction 
kit (Qiagen, Hilden, Germany). Both strands of the PCR products were 
sequenced twice with an ABI Prism 3700 DNA analyzer (Applied Biosys- 
tems, Foster City, CA), using the two PCR primers. The sequences of the 
PCR products were compared with known sequences of the RdRp genes of 
CoVs in the GenBank database. 

Complete genome sequencing. Two complete genomes of porcine 
coronavirus HKU15 (PorCoV HKU15) and one complete genome each of 
white-eye coronavirus HKU16 (WECoV HKU16), sparrow coronavirus 
HKU17 (SpCoV HKU17), magpie robin coronavirus HKU18 (MRCoV 
HKU18), night heron coronavirus HKU19 (NHCoV HKU19), wigeon 
coronavirus HKU20 (WiCoV HKU20), and common moorhen corona- 
virus HKU21 (CMCoV HKU21) were amplified and sequenced using the 
RNA extracted from the original swab specimens as templates. The RNA 
was converted to cDNA by a combined random-priming and oligo(dT)- 
priming strategy. The cDNA was amplified by degenerate primers de- 
signed by multiple alignments of the genomes of other CoVs with com- 
plete genomes available, using strategies described in our previous 
publications (28, 45, 48, 49) and the CoV database CoVDB (20) for se- 
quence retrieval. Additional primers were designed from the results of the 
first and subsequent rounds of sequencing. These primer sequences are 
available on request. The 5’ ends of the viral genomes were confirmed by 
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FIG 1 Phylogenetic analysis of amino acid sequences of the 228-bp fragment 
(excluding primer sequences) of RNA-dependent RNA polymerase (RdRp) of 
CoVs identified from dead wild birds and pigs in the present study. The tree 
was constructed by the neighbor joining method using Kimura correction and 
bootstrap values calculated from 1,000 trees. The scale bar indicates the esti- 
mated number of substitutions per 20 amino acids. The eight genomes 
completely sequenced are shown in bold. PEDV, porcine epidemic diarrhea 
virus (NC_003436); Sc-BatCoV-512, Scotophilus bat coronavirus 512 
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rapid amplification of cDNA ends (RACE) using the 5’/3’ RACE kit 
(Roche, Germany). Sequences were assembled and manually edited to 
produce final sequences of the viral genomes. 

Genome analysis. The nucleotide sequences of the genomes and the 
deduced amino acid sequences of the open reading frames (ORFs) were 
compared to those of other CoVs using EMBOSS needle (http://www.ebi 
.ac.uk). Phylogenetic tree construction was performed using the neighbor 
joining method with ClustalX 1.83. Protein family analysis was performed 
using PFAM and InterProScan (1, 2). Prediction of transmembrane do- 
mains was performed using TMpred and TMHMM (19, 41). 

Estimation of divergence dates. Divergence times for the four genera 
of CoVs, Alphacoronavirus, Betacoronavirus, Gammacoronavirus, and Del- 
tacoronavirus, were calculated using a Bayesian Markov chain Monte 
Carlo (MCMC) approach as implemented in BEAST (Version 1.6.1) as 
described previously (11, 23, 27, 47). One parametric model (Constant 
Size) and one nonparametric model (Bayesian Skyline) tree priors were 
used for the inference. Analyses were performed under the GIR+I+G 
substitution model for RdRp gene sequence data and using both a strict 
and an unrelaxed log-normal-distributed (Ucld) relaxed molecular clock. 
The MCMC run was 5 X 10’ steps long, with sampling every 1,000 steps. 
Convergence was assessed on the basis of the effective sampling size after 
a 10% burn-in using Tracer software version 1.5 (11). The mean time of 
the most recent common ancestor (tMRCA) and the highest posterior 
density regions at 95% (HPD) (i.e., a credible set that contains 95% of the 
sampled values) were calculated, and the best-fitting model was selected 
by a Bayes factor, using marginal likelihoods implemented in Tracer (see 
Table S1 in the supplemental material) (42). Bayesian Skyline under a 
relaxed-clock model with Ucld was adopted for making inferences, as 
Bayes factor analysis indicated that this model fitted the data better than 
other models tested (see Table S1). The trees were summarized in a target 
tree by the Tree Annotator program included in the BEAST package by 


(NC_009657); TGEV, transmissible gastroenteritis virus (NC_002306); FIPV, 
feline infectious peritonitis virus (AY994055); CCoV, canine coronavirus 
(GQ477367); PRCV, porcine respiratory coronavirus (DQ811787); Rh- 
BatCoV-HKU2, Rhinolophus bat coronavirus HKU2 (EF203064); Mi-BatCoV 
1A, Miniopterus bat coronavirus 1A (NC_010437); Mi-BatCoV 1B, Mini- 
opterus bat coronavirus 1B (NC_010436); Mi-BatCoV-HKU8, Miniopterus bat 
coronavirus HKU8 (NC_010438); HCoV-229E, human coronavirus 229E 
(NC_002645); HCoV-NL63, human coronavirus NL63 (NC_005831); HCoV 
OC43, human coronavirus OC43 (NC_005147); BCoV, bovine coronavirus 
(NC_003045); AntelopeCoV, sable antelope CoV (EF424621); GiCoV, giraffe 
coronavirus (EF424622); ECoV, equine coronavirus (NC_010327); PHEV, 
porcine hemagglutinating encephalomyelitis virus (NC_007732); MHV, mu- 
rine hepatitis virus (NC_001846); RCoV, rat coronavirus (NC_012936); 
HCoV-HKU1, human coronaivurs HKU1 (NC_006577); Ty-BatCoV-HKU4, 
Tylonycteris bat coronavirus HKU4 (NC_009019); Pi-BatCoV-HKUS5, Pipist- 
rellus bat coronavirus HKU5 (NC_009020); SARS CoV, SARS-related human 
coronavirus (NC_004718); SARSr-Rh-BatCoV HKU3, SARS-related Rhinolo- 
phus bat coronavirus HKU3 (DQ022305); SARSr CoV CFB, SARS-related 
Chinese ferret badger coronavirus (AY545919); SARSr-CiCoV, SARS-related 
palm civet coronavirus (AY304488); Ro-BatCoV-HKU9, Rousettus bat coro- 
navirus HKU9 (NC_009021); IBV, infectious bronchitis virus (NC_001451); 
IBV-partridge, partridge coronavirus (AY646283); TCoV, turkey corona- 
virus (NC_010800); IBV-peafowl, peafowl coronavirus (AY641576); 
BWCoV-SWI1, beluga whale coronavirus SW1 (NC_010646); ALCCoV, 
Asian leopard cat coronavirus (EF584908); BuCoV HKUL1I, bulbul coro- 
navirus HKU11(FJ376619); ThCoV HKU12, thrush coronavirus HKU12 
(FJ376621); MunCoV HKU13, munia coronavirus HKU13 (FJ376622); 
PorCoV HKU1S5, porcine coronavirus HKU15; WECoV HKU16, white-eye 
coronavirus HKU16; SpCoV HKU17 (TrSp, tree sparrow), sparrow coro- 
navirus HKU17; MRCoV HKU18 (OMR, oriental magpie robin), magpie 
robin coronavirus HKU18; NHCoV HKU19 (BICrNH, black-crowned 
night heron), night heron coronavirus HKU19; WiCoV HKU20 (EuWi, 
Eurasian wigeon), wigeon coronavirus HKU20; CMCoV HKU21, common 
moorhen (CM) coronavirus HKU21. Mu, munia; ChMu, chestnut munia; 
GbTh, gray-backed thrush; ShBu, sooty-headed bulbul; RwBu, red- 
whiskered bulbul; ChBu, chestnut bulbul. 
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TABLE 2 Comparison of genomic features and amino acid identities among CoVs with complete genome sequences available* 


Genome features —_ Pairwise amino acid identity (%) 


WECoV HKU16 SpCoV HKU17 


N 3CLP*° RdRp Hel S N 3CLP*° RdRp Hel S N 


PorCoV HKUI15 
Size G+C 
CoV (bases) content 3CLP7° RdRp Hel S$ 
Alphacoronavirus 
PEDV 28,033 0.42 35.8 48.7 49.3 38.0 
TGEV 28,586 0.38 34.9 49.6 51.6 35.5 
FIPV 29,355 0.38 35.7 49.7 51.2 35.1 
CCoV 29,363 0.38 35.6 49.7 51.6 34.9 
PRCV 27,550 0.37 34.9 49.5 51.6 40.3 
HCoV-229E 27,317 0.38 34.4 49.3 50.6 42.5 
HCoV-NL63 27,553 0.34 35.9 48.8 49.9 38.2 
Rh-BatCoV-HKU2 27,165 0.39 34.4 50.1 51.4 25.0 
Mi-BatCoV 1A 28,326 0.38 33.5 49.0 51.4 35.8 
Mi-BatCoV 1B 28,476 0.39 34.2 48.5 51.1 35.6 
Mi-BatCoV-HKU8 28,773 0.42 33.1 49,3 49.8 35.9 
Sc-BatCoV-512 28,179 0.40 33.8 48.6 49.1 39.0 
Betacoronavirus 
Subgroup A 
HCoV-OC43 30,738 0.37 38.1 51.6 48.3 26.0 
BCoV 31,028 0.37 38.5 51.8 48.4 25.7 
PHEV 30,480 0.37 38.5 51.7 48.3 26.9 
AntelopeCoV 30,995 0.37 38.5 51.8 48.4 25.8 
GiCoV 30,979 0.37 38.8 51.8 48.4 25.9 
ECoV 30,992 0.37 38.5 51,7 49.8 26.0 
MHV 31,357 0.42 38.3 51.9 48.1 26.3 
HCoV-HKU1 29,926 0,32 38.1 51.2 49.3 26.1 
RCoV 31,250 0.41 38.7 51.8 47.9 27.2 
Subgroup B 
SARS CoV 29,751 0.41 34.5 50.7 51.4 26.1 
SARSr-CiCoV 29,728 0.41 34.5 50.7 51.4 26.2 
SARSr-Rh-BatCoV HKU3 29,704 0.41 34.2 50.5 51.4 26.4 
SARSr CoV CFB 29,734 0.41 34.5 50.6 51.4 26.1 
Subgroup C 
Ty-BatCoV-HKU4 30,286 0.38 36.9 51.2 49.8 26.6 
Pi-BatCoV-HKU5 30,488 0.43 35.7 Sel 50.0 26.0 
Subgroup D 
Ro-Bat-CoV HKU9 29,114 0.41 36.4 51.6 51.2 28.4 
Gammacoronavirus 
IBV 27,608 0.38 43.9 54.8 56.6 30.3 
TCoV 27,657 0.38 43.6 54.9 57.1 30.1 
BWCoV-SW1 31,686 0.39 38.8 52.9 52.8 27.1 
Deltacoronavirus 
BuCoV HKU11 26,476 0.39 81.1 88.2 89.4 69.8 
ThCoV HKU12 26,396 0.38 82.1 88.2 89.7 47.9 
MunCoV HKU13 26,552 0.43 82.7 90.1 95.8 71.2 
PorCoV HKUI5 25,421 0.43 
WECoV HKU16 26,027 0.40 76.9 88.1 88.4 61.9 
SpCoV HKU17 26,067 0.45 97.0 97.8 99.2 44.8 
MRCoV HKU18 26,674 0.47 84.3 90.6 96.1 44.4 
NHCoV HKU19 26,064 0.38 54.0 72.5 78.2 41.8 
WiCoV HKU20 26,211 0.39 57.7 71.0 74.9 43.8 
CMCoV HKU21 26,216 0.35 73.6 84.5 84.6 50.1 


choosing the tree with the maximum sum of posterior probabilities (max- 
imum clade credibility) after a 10% burn-in. 

Nucleotide sequence accession numbers. The nucleotide sequences 
of the eight genomes of PorCoV HKU15, WECoV HKU16, SpCoV 
HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and 
CMCoV HKU21 have been lodged within the GenBank sequence data- 
base under accession no. JQ065042 to JQ065049. 


RESULTS 


Animal surveillance and identification of seven novel mamma- 
lian and avian CoVs. A total of 7,140 respiratory and alimentary 
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23.4 37.4 49.3 47.8 38.7 22.4 36.5 48.9 48.9 39.2 24.1 
23.2 34.5 49.4 49.6 36.1 24.5 35.3 49.8 51.4 39.4 23.5 
24.5 35.5 49.6 49.3 36.3 25.1 35.7 49.9 51.1 38.5 25.0 
23.3 35.2 49.4 49.6 35.6 24.0 35.9 49.8 51.4 38.8 23.3 
23.2 34.5 49.3 49.6 40.5 23.2 35.3 49.7 51.4 44.8 23.5 
21.6 35.4 49.0 48.3 42.4 22.5 34.2 49.5 50.2 45.5 23.0 
22.1 38.1 49.2 48.1 40.1 23.0 35.6 49.2 49.6 39.3 22.6 
20.8 34.3 50.0 49.1 25.2 22.3 34.4 50.2 51.1 26.2 20.9 
24.4 35.0 49.4 50.1 35.7 23.3 34.2 49.4 51.1 39.4 25.2 
24.6 35.4 48.8 494 36.1 22.1 34.8 48.8 50.7. 39.1 24.9 
19.4 36.0 49.9 475 36.0 188 33.4 49.6 49.3 38.9 20.4 
24.8 36.0 49.2 47.5 38.7 23.7 34.1 48.7 48.8 41.3 25.2 


22.2 38.9 Bilis} 48.6 25.9 23.2 37.8 51.8 48.3 26.9 22.4 
22.9 38.8 51.7 48.6 25.8 21.7 38.5 51.8 48.4 26.7 22.8 
22.1 38.1 51.6 48.6 26.1 23.1 38.5 51.6 48.3.9 27.2 22.3 
22.9 38.8 51.7 48.5 25.6 21.7 38.5 51.8 48.4 27.0 22.1 
22.9 38.8 51.7 48.5 25.7 21.7 38.8 51.8 48.4 27.2 22.1 
23.9 38.8 51.6 49.0 264 22.6 38.5 51.7 49.9 26.5 24.0 
24.3 39.0 51.3 48.5 26.1 24.0 38.3 51.8 48.3 26.5 25.3 
25.2 38.0 51.4 48.2 264 24.8 37.9 51.3 49.4 25.7 26.0 
24.5 39.5 51.4 48.3 27.0 24.3 38.5 51.7 48.1 25.5 25.1 


26.5 36.1 50.3 50.6 27.9 24.7 34.2 51.1 51.6 25.3 25.6 
26.5 36.1 50.3 50.6 28.0 24.7 34.2 51.1 51.6 25.2 25.6 
25.2 35.8 50.3 50.8 26.2 24.3 33.9 o11 51.6 25.6 24.9 
265 36.1 50.2 50.6 28.0 24.7 34.2 51.0 51.6 25.5 25.6 


25.1 36.6 51.0 49.4 26.1 24.7 36.9 51.5 49.7 27.0 26.2 
25.6 37.8 50.3 49.0 25.5 25.7 35.4 51.4 49.8 27.2 25.3 


25.1, 39.2 52.6 50.1 26.5 23.1 36.4 51.7 50.9 26.6 24.9 


30.0 42.6 54.6 54.5 29.9 30.8 44.2 54.9 56.6 27.6 28.9 
29.2 43.3 54.5 55.3 30.3 29.6 43.9 55.0 57.1 29.5 29.5 
32.1 39.5 52.9 51.6 28.3 31.1 388 52.9 52.8 28.5 31.9 


74.8 82.4 90.9 96.0 62.5 73.2 80.8 88.2 89.6 43.5 75.1 
79.7 83.1 89.5 94.7 47.8 81.0 81.8 88.2 89.9 46.7 79.4 
76.8 76.5 87.9 89.1 61.3 74.6 83.4 90.1 96.0 43.8 78.8 
76.9 88.1 88.4 61.9 75.8 97.0 97.8 99.2 448 96.8 
75.8 77.2 88.3 88.6 464 76.4 
96.8 77.2 88.3 88.6 464 76.4 
77.9 77.2 87.3 89.1 45.5 75.4 84.9 91.0 96.3 68.1 79.0 
52.2 55.0 71.7 76.6 42.1 49.6 52.8 72.3 78.4 47.2 51.9 
52.4 59.6 71.5 74.3, 43.4 50.6 58.0 71.0 75.1 45.8 53.2 
62.0 76.9 84.8 90.5 51.5 64.3 73.0 84.7 84.9 46.0 63.5 


(Continued on following page) 


specimens from 3,298 dead wild birds, 221 chickens, and 3,137 
mammals were obtained (Table 1). RT-PCR for a 440-bp frag- 
ment in the RdRp genes of CoVs was positive in specimens 
from 17 pigs and 35 dead wild birds. Sequencing results sug- 
gested the presence of seven novel CoVs (Fig. 1 and Table 1). 
These seven novel CoVs were most closely related to our re- 
cently described BuCoV HKU11, ThCoV HKU12, and Mun- 
CoV HKU13, sharing <66% nucleotide identity with all other 
known CoVs (Fig. 1). No positive results were obtained from 
any of the 15 Asian leopard cats, 434 bats, 230 cats, 47 cattle, 
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TABLE 2 (Continued) 


Pairwise amino acid identity (%) 


MRCoV HKUI18 NHCoV HKU19 WiCoV HKU20 CMCoV HKU21 
3CLP*° RdRp _—-Hel S N 3CLP*° RdRp Hel S N 3CLP*° = RdRp _ Hel S N 3CLP*° =RdRp_—s Hel S N 
36.8 48.9 49.1 40.8 22.3 37.2 48.8 47.9 36.3 20.6 39.7 50.0 48.5 37.9 21.1 38.1 49.6 46.7 39.4 21.7 
34.6 49.6 50:7. 39.9 23.7 32.9 50.8 50.8 37.1 22.3 37.3 50.2 50.2 36.7 23.7 33.7 49.3 494 364 22.9 
35.5 50.1 50.4 38.9 22.2 32.4 51.2 50.7. 36.9 20.8 37.9 50.1 50.0 365 22.5 34.5 49.7 49.1 36.0 22.1 
35.3 49.6 50.7. 39.3. 23.0 = 32.9 50.8 50.7. 36.4 23.6 37.3 49.9 50.7. 36.3. 23.7 33.3 49.1 494 35.2 23.4 
34.6 49.5 50.7. 444 228 32.9 50.7 50.8 414 22.0 37.3 50.3 50.2. 41.9 23.7 33.7 49.1 49.4 40.1 21.8 
34.8 49,3 498 440 22.5 34.6 49.2 490 39.9 204 36.3 49.6 496 44,1 215 34.8 48.6 47.5 43.2 22.8 
36.9 49,2 495 39.3 246 34.8 48.9 490 362 21.5 39.1 50.6 48.6 388 23.5 36.9 49.5 47.2 38.8 20.6 
34.4 49.8 50.3. 25.9 21.7 34.3 50.5 49.7 25.1 20:9. (35.5 50.6 49.1 26.3 25.2 34.1 50.1 49.1 27.3, 22.5 
34.2 49.0 51.4 382 245 32.6 50.5 48.7 35.7 23.8 34.7 49.9 489 364 22.55 33.4 47.9 49.1 38.6 22.4 
33.5 48.8 51.1 38,2 23.8. 31.9 49.8 48.2 35.7 22.3 35.7 49.6 48.9 35.7 24.1 32.8 47.7 48.3 38.2 22.8 
34.4 49,3 49.3 40.4 19.8 34.6 50.1 48.1 37.2 22.6 36.7 50.5 48.8 37.6 22.3 36.0 48.5 474 37.0 21.6 
35.4 48.3 48.8 41.1 23.4 34.9 49.0 478 368 22.3 37.3 48.9 49.1 38.1 23.1 34.7 48.7 478 39.6 24.3 
37.5 51.3 48.3 264 21.0 34.1 54.5 48.4 25.9 22.5 38.7 51.8 48.4 25.7 204 37.8 51.5 49.0 25.6 24.2 
37.8 51.5 484 269 240 34.4 54.5 48.5 25.6 23.4 38.3 51.7 48.5 260 21.5 38.2 51.6 49.0 25.7 23.4 
37.8 51.4 48.3 269 22.2 34.4 54.5 48.5 26.1 22.7 38.7 51.7 48.5 27.1 216 38.2 51.6 490 25.4 24.1 
37.8 51.4 484 27.0 240 34.4 54.4 48.5 26.1 23.9 38.3 51.7 485 262 21.5 38.2 51.6 49.2 25.8 23.4 
37.8 51.4 484 27.0 240 34.4 54.4 48.5 25.7 23.9 38.3 51.7 48.5 265 21.5 38.5 51.6 49.2 25.9 23.4 
37.8 51.4 498 275 23.5 34.4 54.6 48.5 25.3 246 38.3 oS 48.5 269 22.0 38.2 51.6 49.7 25.6 24.9 
37.6 51.9 48.3 263 24.2 35.0 53.6 475 25.3 246 39.6 50.8 479 27.1 24.0 39.2 51.2 48.6 26.0 24.6 
36.4 51.4 488 264 26.0 36.3 54.4 474 254 24.7 38.1 50.9 48.5 25.8 22.7 38.3 51.2 478 25.0 25.4 
38.2 51.8 48.2 25.8 25.0 35.0 53.6 474 243 25.2 39.9 50.7 47.7 274 24.1 38.3 51.0 48.4 264 23.5 
34.2 50.8 514 254 26.2 32.1 50.5 50.2 26.3 22.7 34.8 49.8 50.3 26.9 24.3 32.9 50.8 51.0 27.3 24.8 
34.2 50.8 514 25.5 26.2 32.1 50.5 50.2 26.2 22.7 34.8 49.8 50.3. 27.0 24.3 32.9 50.8 51.0 27.1 24.8 
33.9 50.9 514 260 25.7 32.1 50.4 50.6 25.6 23.0 34.8 49.7 50.5 26.0 23.5 32.6 50.6 51.30 27.2 24.1 
34.2 50.7 514 25.2 26.2 32.1 50.4 50.2 25.9 22.7 34.8 49.9 50.3. 26.8 24.3 32.9 50.7 51.0 27.2 24.8 
35.7 51.0 49.7 27.3 25.9 32.7 51.3 499 27.3 244 35.8 50.9 499 264 244 35.6 51.9 489 268 25.1 
35.0 51.1 49.7 26.3 26.2 33.7 50.9 496 26.2 25.6 34.6 51.2 498 25.3 24.7 36.0 50.9 49.0 25.7 26.1 
35.8 51.9 50.9 27.7 25.0 33.8 52.2 50.9 27.0 22.8 35.0 51.4 492 27.2 22.9 36.9 52.3 514 27.8 23.3 
43.3 54.3 56.2 284 30.3 43.6 53.6 54.8 28.7 29.7 47.1 52.4 55.4 294 29.2 41.7 53.9 55.1 29.4 28.1 
42.3 54.4 57.1 30.3. 29.2 44.3 53.8 55.3 30.3 30.0 46.2 52.8 55.4 30.3 29.9 41.3 53.6 55.8 30.6 28.2 
37.2 52.3 52.2 27.33 316 41.1 52.8 54.2 274 30.2 42.0 52.2 51.1 27.2 30.0 38.8 52.1 52.7 28,2. 31.5 
79.5 88.3 90.4 445 71.9 57.0 72.3 76.8 41.1 50.6 58.3 70.8 75.4 43.3 514 77.5 84.8 91.0 51.8 60.7 
81.4 86.8 89.9 45.8 76.7 57.3 71.9 764 43.6 49.4 57.7 71.3 74.5 43.6 49.6 78.2 84.4 90.5 46.2 63.3 
94.5 94.6 98.0 46.1 87.5 53.1 72.9 78.0 414 53.4 55.4 71.7 75.7 44.0 53.2 72.0 84.7 85.4 52.2 64.4 
84.3 90.6 96.1 444 77.9 54.0 72.5 78.2 41.8 52.2 57.7 71.0 74.9 43.8 524 73.6 84.5 84.6 50.1 62.0 
772 87.3 89.1 45.5 75.4 55.0 77 76.6 42.1 49.6 59.6 71.5 74.3 434 50.6 76.9 84.8 90.5 51.5 64.3 
84.9 91.0 96.3 68.1 79.0 52.8 72.3 78.4 47.2 51.9 58.0 71.0 75.1 45.8 53.2 73.0 84.7 84.9 46.0 63.5 
54.0 72.5 77.7 464 53.8 56.4 71.2 75.1 46.3 53.1 73.3 85.1 84.8 45.7 63.9 
54.0 72.5 77.7 464 53.8 58.3 69.3 75.4 41.0 545 55.5 71.9 77.6 43.6 54.5 
56.4 71.2 75.1 46.3 53.1 58.3 69.3 75.4 41.0 54.5 58.3 70.8 76.4 44.1 57.0 
73.3 85.1 84.8 45.7 63.9 55.5 71.9 77.6 43.6 54.5 58.3 70.8 76.4 44.1 57.0 


* Comparison of genomic features of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 and other 
CoVs with complete genome sequences available and of amino acid identities between the predicted 3CL’*°, RNA-dependent RNA (RdRp), helicase (Hel), S, and N proteins of 
PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21 and the corresponding proteins of other CoVs. 
PEDV, porcine epidemic diarrhea virus; TGEV, porcine transmissible gastroenteritis virus; FIPV, feline infectious peritonitis virus; CCoV, canine coronavirus; PRCV, porcine 
respiratory coronavirus; HCoV-229E, human coronavirus 229E; HCoV-NL63, human coronavirus NL63; Rh-BatCoV-HKU2, Rhinolophus bat coronavirus HKU2; Mi-BatCoV 1A, 
Miniopterus bat coronavirus 1A; Mi-BatCoV 1B, Miniopterus bat coronavirus 1B; Mi-BatCoV-HKU8, Miniopterus bat coronavirus HKU8; Sc-BatCoV-512, Scotophilus bat 
coronavirus 512; HCoV OC43, human coronavirus OC43; BCoV, bovine coronavirus; PHEV, porcine hemagglutinating encephalomyelitis virus; AntelopeCoV, sable antelope 
coronavirus; GiCoV, giraffe coronavirus; ECoV, equine coronavirus; MHV, murine hepatitis virus; HCoV-HKU1, human coronavirus HKU1; RCoV, rat coronavirus; SARS CoV, 
SARS-related human coronavirus; SARSr-CiCoV, SARS-related palm civet coronavirus; SARSr-Rh-BatCoV HKU3, SARS-related Rhinolophus bat coronavirus HKU3; SARSr CoV 
CEB, SARS-related Chinese ferret badger coronavirus; Ty-BatCoV-HKU4, Tylonycteris bat coronavirus HKU4; Pi-BatCoV-HKU5, Pipistrellus bat coronavirus HKU5; 
Ro-BatCoV-HKU9, Rousettus bat coronavirus HKU9; IBV, infectious bronchitis virus; TCoV, turkey coronavirus; BWCoV-SW1, Beluga whale coronavirus SW1; BuCoV HKU11, 
bulbul coronavirus HKU11; ThCoV HKU12, thrush coronavirus HKU12; MunCoV HKU13, munia coronavirus HKU13. 
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FIG 2 Genome organization of members in Deltacoronavirus. ORFs downstream of S gene are magnified to show the differences among the genomes of the 10 
CoVs. Papain-like protease (PL?"°), chymotrypsin-like protease (3CL?”°), and RNA-dependent RNA polymerase (RdRp) are represented by orange boxes. Spike 
(S), envelope (E), membrane (M), and nucleocapsid (N) are represented by green boxes. Putative accessory proteins are represented by blue boxes. The seven 


CoVs discovered in this study are shown in bold. 


221 chickens, 231 dogs, 1,387 humans, 235 monkeys, and 389 
rodents tested (Table 1). 

Genome organization and coding potential of the seven 
novel mammalian and avian CoVs. Complete genome sequence 
data of two strains of PorCoV HKU15 and one complete ge- 
nome each of WECoV HKU16, SpCoV HKU17, MRCoV 
HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV 
HKU21 were obtained by assembly of the sequences of the 
RT-PCR products from the RNA extracted from the corre- 
sponding individual specimens. 

The size of the genomes of the seven novel CoVs ranged from 
25,416 bases (PorCoV HKU15) to 26,674 (MRCoV HKU18) and 
their G+C contents ranged from 35% (CMCoV HKU21) to 47% 
(MRCoV HKU18) (Table 2). Their genome organizations are 
similar to those of other CoVs, with the characteristic gene order 
5'-replicase ORF lab, spike (S), envelope (E), membrane (M), nu- 
cleocapsid (N)-3' (Fig. 2 and Table 3). Both 5’ and 3’ ends contain 
short untranslated regions. The replicase ORF lab occupies 18.620 
to 18.887 kb of the genomes (Table 3). This ORF encodes a num- 
ber of putative proteins, including nsp3 [which contains the pu- 
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tative papain-like protease (PL?™°)], nsp5 [putative chymo- 
trypsin-like protease (3CLP*°)], nsp12 (putative RdRp), nsp13 
(putative helicase), and other proteins of unknown functions. No- 
tably, the amino acids upstream to the putative cleavage sites at 
nsp2/nsp3, nsp3/nsp4, and nsp4/nsp5 are all AG, AG, and LQ for 
PorCoV HKUI5, WECoV HKU16, SpCoV HKU17, MRCoV 
HKU18, and CMCoV HKU21; however, those at nsp2/nsp3 are 
VG and DG, those at nsp3/nsp4 are TG and GG, and those at 
nsp4/nsp5 are VQ for NHCoV HKU19 and WiCoV HKU20 (see 
Table S3 in the supplemental material). 

The seven novel CoVs display similar genome organizations 
and differ only in the number of ORFs downstream of N (Fig. 2). 
Their transcription regulatory sequences (TRSs) conform to the 
consensus motif 5’-ACACCA-3’ (Table 3), which appears to be 
unique to members of the genus Deltacoronavirus. Interestingly, 
similar to BuCoV HKU11, ThCoV HKU12, and MunCoV 
HKU13, the perfect TRSs of S in the genomes of the seven novel 
CoVs were separated from the corresponding AUG by 80 to 145 
bases (Table 3). This is in contrast to the relatively small number of 
bases between the TRSs for S and the corresponding AUG (range: 
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TABLE 3 Coding potential and putative transcription regulatory sequences of CoV genomes* 


Putative TRS 
CoV ORF Location (nt) Length (nt) Length (aa) Frame TRS location (nt) TRS sequence(s) (distance in bases to AUG)? 
PorCoV HKU15 lab 540-19342 18,803 6,268 +342 75 ACACCA(459)AUG 
S 19324-22806 3,483 1,161 “EL 19178 ACACCA(145)AUG 
E 22800-23051 252 84 +3 22777 ACACCG(17)AUG 
M 23044-23697 654 218 Sal 23018 ACACCA(20)AUG 
NS6 23697-23981 285 95 +3 23645 ACACCA(46)AUG 
N 24002-25030 1,029 343 +2 23989 ACACCA(7)AUG 
NS7 24096-24698 603 201 +3 24008 GCACCA(82)AUG 
WECoV HKU16 lab 511-19397 18,887 6,296 +#1,;4+3 66 ACACCA(439)AUG 
S 19379-22918 3,540 1,180 +2 19233 ACACCA(140)AUG 
E 22912-23160 249 83 +1 22886 ACACCA(20)AUG 
M 23153-23809 657 219 +2 23130 ACACCA(17)AUG 
NS6 23809-24090 282 94 ae 23768 ACAUCA(35)AUG 
N 24115-25158 1,044 348 ol 24101 ACACCA(8)AUG 
NS7a 24143-24811 669 223 +2, 24101 ACACCA(36)AUG 
NS7b 25139-25270 132 44 +2 25039 AAACCA(94)AUG 
SpCoV HKU17 lab 520-19352 18,833 6,278 #1, +3 57 ACACCA(452)AUG 
S 19334-22954 3,621 1,207 +2 19188 ACACCA(140)AUG 
E 22948-23196 249 83 +E 22925 ACACCG(17)AUG 
M 23189-23842 654 218 +2 23166 ACACCA(17)AUG 
NS6 23842-24129 288 96 +1 23790 ACACCA(46)AUG 
N 24150-25178 1,029 343 +3 24137 ACACCA(7)AUG 
NS7a 25189-25623 435 145 al 25179 ACACCA(4)AUG 
NS7b 25539-25751 213 71 +3 25523 ACUCCA(10)AUG 
MRCoV HKU18 lab 596-19356 18,761 6,254 +2, +1 64 ACACCA(526)AUG 
Ss 19338-22991 3,654 1,218 +3 19192 ACACCA(140)AUG 
E 22985-23233 249 83 +2 22945 ACACCG(34)AUG 
M 23226-23882 657 219 +3 23203 ACACCA(17)AUG 
NS6 23882-24172 291 97 +2 23857 ACGCCA(19)AUG 
N 24355-25395 1,041 347 el 24340 ACACCA(9)AUG 
NS7a 25407-25580 174 58 +3 25396 ACACCA(5)AUG 
NS7b 25561-25932 372 124 +1 
NS7c 25941-26195 255 85 +3. 25910 ACACCA(25)AUG 
NHCoV HKU19 lab 482-19323 18,842 6,281 +2,4+1 67 ACACCG(409)AUG 
S 19305-23069 3,765 1,255 +3 19156 ACACCG(143)AUG 
E 23069-23317 249 83 +2 23013 ACACCA(50)AUG 
M 23310-23960 651 217 +3 23211 ACACCG(93)AUG 
NS6 23960-24238 279 93 +2 23951 ACACCU(3)AUG 
N 24248-25276 1,029 343 +2 24231 ACACCU(8)AUG 
NS7a 25277-25573 297 99 +2 25248 ACACCG(23)AUG 
NS7b 25583-25876 294 98 +2 25560 ACACCA(17)AUG 
WiCoV HKU20 lab 219-18838 18,620 6,207 #3, 2. 60 ACACCA(153)AUG 
S 18817-22455 3,639 1,213 Boal 18731 ACACCU(80)AUG 
E 22455-22715 261 87 +3 22380 ACACCA(69)AUG 
M 22708-23358 651 217 al 22597 ACACCG(105)AUG 
NS6 23358-23630 273 91 +3 
N 23646-24698 1,053 351 +3 23631 ACACCA(9)AUG 
NS7a 24695-24928 234 78 +2 24609 AAACCA(80)AUG 
NS7b 25218-25466 249 83 +3 25177 ACACCG(35)AUG 
NS7c 25450-25716 267 89 + 25444 ACACCGAUG 
NS7d 25752-25952 201 67 +3 25735 AAACCU(11)AUG 
CMCoV HKU21 lab 478-19103 18,626 6,209 #13) 63 ACACCA(409)AUG 
S 19085-22729 3,645 1,215 +2 18939 ACACCA(140)AUG 
E 22723-22971 249 83 =e 22697 ACACCA(20)AUG 
M 22973-23779 807 269 +2 22938 ACACCA(29)AUG 
NS6 23779-24024 246 82 +e 23727 ACACCA(46)AUG 
N 24052-25107 1,056 352 + 24039 ACACCG(7)AUG 
NS7a 25107-25379 273 91 +3 25036 ACACCU(65)AUG 
NS7b 25391-25576 186 62 +2 25379 ACACCU(6)AUG 
NS7c 25500-25916 A417 139 +2 


* PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21. aa, amino acid; nt, nucleotide. 
» Boldface indicates putative TRS sequences. The nucleotide variations are in italic. 
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CoV Acc. No. 3°’-UTR stem loop structure 

IBV NC_001451 27,471 CAGTGCCGGGGCCACGCGGAGTACGATCGAGGGTACAGCACTA 27,513 
SARS CoV NC_004718 29,584 TTTCATCGAGGCCACGCGGAGTACGATCGAGGGTACAGTGAAT 29,626 
SARSr-Rh-BatCoV HKU3 DQ022305 29,561 TTTCACCGAGGCCACGCGGAGTACGATCGAGGGTACAGTGAAT 29,603 
BuCoV HKUI1 FJ376619 26,267 ATGTGCCGAGGCCACGCGGAGTACGATCGAGGGTACAGCACAA 26,309 
ThCoV HKU12 FJ376621 26,173 ATATGCCGAGGCCACGCGGAGTACGATCGAGGGTACAGCATAA 26,215 
MunCoV HKU13 FJ376622 26,329 ATGTGTCGAGGCCACGCGGAGTACGATCGAGGGTACAGCACAA 26,371 
ALCCoV EF584908 12,603 ATATGCCGAGGCCACGCGGAGTACGATCGAGGGTACAGCATAA 12,645 
PorCoV HKU15 JQ065042 25,214 ATATGCCGAGGCCACGCGGAGTACGATCGAGGGTACAGCATAA 25,256 
WECoV HKU16 JQ065044 25,819 TTGCACCGAGGCCACGCGGAGTACGATCGAGGGTACAGTGCAC 25,861 
SpCoV HKU17 JQ065045 25,859 ATATGCCGAGGCCACGCGGAGTACGATCGAGGGTACAGCATAA 25,901 
MRCoV HKU18 JQ065046 26,466 ATGTGCCGAGGCCACGCGGAGTACGATCGAGGGTACAGCACAA 26,508 
CMCoV HKU21 JQ065049 26,007 ATGAACCGAGGCCACGCGGAGTACGATCGAGGGTACAGTTCAA 26,049 


TER RRR OR OR Kk 


FIG3 Multiple alignments of conserved s2m of infectious bronchitis virus (IBV), SARS-related human coronavirus (SARS CoV), SARS-related Rhinolophus bat 
coronavirus HKU3 (SARSr-Rh-BatCoV HKU3), BuCoV HKU11, ThCoV HKU12, MunCoV HKU13, Asian leopard cat coronavirus (ALCCoV), PorCoV 
HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, and CMCoV HKU21. Identical nucleotides are marked by asterisks. Acc. No., accession no. 


from 0 bases in HCoV-NL63, Rhinolophus bat coronavirus 
HKU2 [Rh-BatCoV-HKU2], HCoV-HKU1, bovine coronavirus 
[BCoV], HCoV-OC43, mouse hepatitis virus [MHV], porcine 
hemagglutinating encephalomyelitis virus, SARS-CoV, and 
SARS-related Rhinolophus bat coronavirus HKU3 [SARSr-Rh- 
batCoV HKU3] to 52 bases in infectious bronchitis virus [IBV]) in 
members of Alphacoronavirus, Betacoronavirus, and Gammacoro- 
navirus. Similar to BuCoV HKU11, ThCoV HKU12, and 
MunCoV HKU13, the genomes of the seven novel CoVs have 
putative PL?’®, which are homologous to PL2?"° of Alphacorona- 
virus and Betacoronavirus subgroup A and PL?" of Betacoronavi- 
rus subgroups B, C, and D and Gammacoronavirus (Fig. 2). Similar 
to BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13, one 
ORF (NS6) is found between M and N of the genomes of the seven 
novel CoVs. On the other hand, one ORF (NS7) is present over- 
lapping with N in PorCoV HKU15, two ORFs (NS7a and 7b) are 
present overlapping or downstream of N in WECoV HKU16, 
SpCoV HKU17, and NHCoV HKU19, three ORFs (NS7a, 7b, and 
7c) are present downstream of N in MRCoV HKU18 and CMCoV 
HKU21, and four ORFs (NS7a, 7b, 7c, and 7d) are present over- 
lapping or downstream of N in WiCoV HKU20. For NS7 of 
PorCoV, the presence of an imperfect TRS (GCACCA) and its 
relatively high Ka/K, ratio (number of nonsynonymous substitu- 
tions per nonsynonymous site/number of synonymous substitu- 
tions per synonymous site) of 1.046 (data not shown) implied that 
this ORF may not be expressed. BLAST search revealed no amino 
acid similarities between these putative nonstructural proteins 
and other known proteins, and no functional domain was 
identified by PFAM and InterProScan, except that NS7a of 
NHCoV HKU19 was found to be homologous to the NS7a of 
BuCoV HKU11, ThCoV HKU12, and MunCoV HKU13. NS7b 
of WiCoV HKU20 and CMCoV HKU21, and NS7d of WiCoV 
HKU20, were also found to be homologous to the NS3b of IBV 
and hypothetical protein of goose coronavirus, respectively. 
Transmembrane helices, predicted by TMHMM and TMpred, in 
putative accessory proteins downstream to the N genes in the 
genomes of SpCoV HKU17, MRCoV HKU18, NHCoV HKU19, 
WiCoV HKU20, and CMCoV HKU21 are listed in Table $4 in the 
supplemental material. Each of the genomes of PorcoV HKU15, 
WECoV HKU16, SpCoV HKU17, MRCoV HKU18, and CMCoV 
HKU21 contains a stem-loop II motif (s2m) (residues 25,220 to 
25,251, 25,825 to 25,856, 25,865 to 25,896, 26,472 to 26,503, and 
26,013 to 26,044, respectively), a conserved RNA element down- 
stream of N and upstream of the poly(A) tail, similar to those in 
IBV, TCoV, SARSr-Rh-BatCoV, and SARS-CoV, as well as other 
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CoVs discovered in Asian leopard cat, graylag geese, feral pigeons, 
and mallards, for which complete genomes are not available (Fig. 
3) (14, 21, 38). 

Comparison of the amino acid identities of the seven con- 
served replicase domains for species demarcation (ADRP, nsp5 
[3CLP*°], nsp12 [RdRp], nsp13 [Hel], nsp14 [ExoN], nsp15 
[NendoU], and nsp16 [O-MT]) (8) among the 10 deltacoronavi- 
ruses is shown in Table S5 in the supplemental material. In all the 
seven domains, the amino acid sequences of PorCcoV HKU15 and 
SpCoV HKU17 showed more than 90% identity, indicating that 
these two coronaviruses should be subspecies of the same species. 

Phylogenetic analyses. The phylogenetic trees constructed us- 
ing the nucleotide sequences of the 3CL?"°, RdRp, Hel, S, and N of 
the seven novel CoVs and other CoVs are shown in Fig. 4 and the 
corresponding pairwise amino acid identities are shown in Table 
2. For all five genes, the seven novel CoVs possessed higher amino 
acid identities to each other and BuCoV HKU11, ThCoV HKU12, 
and MunCoV HKU13 than to any other known CoVs with complete 
genomes available (Table 2). In all five trees, the seven novel CoVs 
were clustered with BuCoV HKU11, ThCoV HKU12, and MunCoV 
HKU13 (Fig. 4). For Hel, S, and N, PorCoVs were also clustered with 
a CoV found in Asian leopard cat (10), for which the sequences of 
these genes were available (Fig. 4). There were <2% base differences 
between the Hel, S, and N genes of PorCoV and those of the Asian 
leopard cat coronavirus. Based on both phylogenetic tree analyses 
and amino acid differences, the seven novel CoVs as well as BuCoV 
HKU11, ThCoV HKU12, and MunCoV HKU13 should belong to 
the same genus, Deltacoronavirus. 

Estimation of divergence dates. Using the Bayesian Skyline 
under a relaxed-clock model with an uncorrelated log-normal dis- 
tribution, the mean evolutionary rate of CoVs was estimated at 
1.3 X 10° * nucleotide substitutions per site per year for the RdRp 
gene. Molecular clock analysis using the RdRp gene showed that 
the tMRCA ofall CoVs was estimated at ~8100 BC (HPDs, 20607 
to 974 BC), that of Alphacoronavirus at ~2400 BC (HPDs, 7659 to 
722 BC), that of Betacoronavirus at ~3300 BC (HPDs, 9713 to 447 
BC), that of Gammacoronavirus at ~2800 BC (HPDs, 8840 to 700 
BC), and that of Deltacoronavirus at ~3000 BC (HPDs, 9073 to 
555 BC) (Fig. 5). 


DISCUSSION 

The diversity of CoVs in birds is comparable to that observed in 
bats. In the last 7 years, we and others have demonstrated a 
previously unrecognized diversity of CoVs in bats (4, 6, 23, 25, 
26, 28, 43). More than 10 CoVs were discovered in bats, with at 


Journal of Virology 


ysanB Aq GLog ‘8 yore UO /Bio‘wse'IAl//:dy}y Wood papeojumoqg 


Clr” 


1000 -Mi-BatCoV 1A 


Alphacoronavirus 
HCoV-229E 


RhBatCoV-HKU2 
FIPV 
CCoV 


748'TGEV 

SARSr-CiCoV 

ISARS CoV 

ISARSr CoV CFB 
SARSr-Rh-BatCoV HKU3 
Ro-BatCoV-HKUS 


Betacoronavirus 


Gammacoronavirus 


IBV-partridge 
-WiCoV HKU20-9243 EuWi/Hong Kong/2008 
NHCoV HKU19-6918 BICrNH/Hong Kong/2007 
CMCoV HKU21-8295 CM/Hong Kong/2007 
'WECoV HKU16-6847 WE/Hong Kong/2007 
BuCoV HKU11 
ThCoVHKU12 
MRCoV HKU18-chu3 OMR/Hong Kong/2007 
MunCoV HKU13 
SpCoV HKU17-6124 TrSp/Hong Kong/2007 
1000 PorCoV HKU15-44/Hong Kong/2009 

1000 !PorCoV HKU15-155/Hong Kong/2010 


885 


Deltacoronavirus 


mM 
ol 


Novel Mammalian and Avian Deltacoronavirus 


RdRp 


1000 - Mi-BatCoV 1B 
Mi-BatCoV 1A 
Mi-BatCoV-HKU8 
RhBatCoV-HKU2 
HCoV-NL63 
HCoV-229E 


Alphacoronavirus 


659 


‘SARSr-Rh-BatCoV HKU3 

Ro-BatCoV-HKUS 

Pi-BatCoV-HKUS 

Betacoronavirus 


ss3!AntelopeCoV 
BWCoV-SW1 
IBV-partridge 


Gammacoronavirus 


856 LIBV-peafowl 
NHCoV HKU19-6918 BICrNH/Hong Kong/2007 
WiCoV HKU20-9243 EuWi/Hong Kong/2008 
-CMCoV HKU21-8295 CM/Hong Kong/2007 
BuCoVHKU11 
ThCoV HKU12 
-WECoV HKU16-6847 WE/Hong Kong/2007 
MRCoV HKU18-chu3 OMR/Hong Kong/2007 
MunCoV HKU13, 
9x4 SpCoV HKU17-6124 TrSp/Hong Kong/2007 
009 | PorCoV HKU15-44/Hong Kong/2009 

1000 'PorCoV HKU15-155/Hong Kong/2010 


Deltacoronavirus 
1000 


342 


— 
00s 


FIG 4 Phylogenetic analyses of 3CL?"°, RdRp, helicase (Hel), S, and N proteins of PorCoV HKU15, WECoV HKU16, SpCoV HKU17, MRCoV HKU18, 
NHCoV HKU19, WiCoV HKU20, and CMCoV HKU21. The trees were constructed by using the neighbor joining method using Kimura correction and 
bootstrap values calculated from 1,000 trees. Two hundred ninety-five, 892, 590, 802, and 249 amino acid positions in 3CL?*°, RdRp, Hel, S, and N, 
respectively, were included in the analyses. The trees were midpoint rooted. For 3CL?"° and S, the scale bar indicates the estimated number of substitutions 
per 10 amino acids. For RdRp and Hel, the scale bar indicates the estimated number of substitutions per 20 amino acids. For N, the scale bar indicates the 
estimated number of substitutions per 5 amino acids. Viruses characterized in this study are in bold. Virus name abbreviations are the same as those in 


the Fig. 1 legend. 


least nine present in our locality, and complete genome se- 
quences are available for eight, which includes SARSr-Rh- 
BatCoV HKU3, Rh-BatCoV-HKU2, Miniopterus bat coronavi- 
rus 1, Miniopterus bat coronavirus HKU8, Scotophilus bat 
coronavirus 512, Tylonycteris bat coronavirus HKU4, Pipistrel- 
lus bat coronavirus HKU5, and Rousettus bat coronavirus 
HKU9 (4, 6, 25, 26, 43). Due to the similarities between bats 
and birds, such as their abilities to fly and high species diversity, 
we hypothesized that there should be previously unrecognized 
CoVs in birds. In our previous study and the present one, we 
demonstrated that there are at least nine CoVs, in addition to 
IBV and its close relatives, in birds (49). Potentially novel CoVs 
in Gammacoronavirus were also observed in another study, al- 
though complete genome sequences are not available and 
therefore detailed genomic and phylogenetic analysis are not 
possible (35). The nine CoVs discovered in the present and 
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previous studies were found in birds of nine different families, 
showing host specificity. This phenomenon of host specificity 
is similar to that observed in bats, in which different genera are 
hosts of different CoVs (26, 45, 51, 52). We speculate that this 
diversity and host specificity of bat and bird CoVs is due to the 
large variety of species in bats and birds, giving rise to a large 
variety of cell types and receptors for the different CoVs to 
attach and replicate. 

The presence of a huge diversity of bat CoVs in Alphacorona- 
virus and Betacoronavirus but not Gammacoronavirus and 
Deltacoronavirus and a huge diversity of bird CoVs in Gamma- 
coronavirus and Deltacoronavirus but not Alphacoronavirus and 
Betacoronavirus supports our model of CoV evolution, in which 
bats are the gene source of Alphacoronavirus and Betacoronavirus 
and birds the gene source of Gammacoronavirus and Deltacorona- 
virus (Fig. 6) (52). It is not known whether the first CoVs occurred 
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FIG 4 continued 


in bats and jumped to birds or vice versa. In the bat CoV lineage, 
the bat CoV jumped to another species of bat, giving rise to Al- 
phacoronavirus and Betacoronavirus. These bat CoVs in turn 
jumped to other bat species and other mammals, including hu- 
mans, with each interspecies jumping evolving dichotomously. As 
for the bird CoV lineage, the bird CoV jumped to another species 
of bird, giving rise to Gammacoronavirus and Deltacoronavirus. 
These bird CoVs in turn jumped to other bird species and occa- 
sionally to some mammalian species, such as whale and pig, with 
each interspecies jumping evolving dichotomously. Although 
PorCoV HKU15 was closely related to a CoV previously found in 
Asian leopard cats and Chinese ferret badgers, further experi- 
ments are warranted to confirm whether these viruses really rep- 
licate in the corresponding animals. Of note is that the estimation 
of divergence time was based on a relaxed-clock assumption with 
no recombination among the genomes. Since CoVs have a ten- 
dency to recombine, the estimated divergence time gives only a 
rough approximation of the actual divergence time. When more 
complete genomes of CoVs in the four different genera at different 
time points are available, such divergence time estimation can be 
performed using multiple gene loci to achieve more accurate esti- 
mation. 
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Both avian and mammalian CoVs are members of Delta- 
coronavirus, with similar genome characteristics and struc- 
tures. In all the 10 members of Deltacoronavirus with complete 
genome sequences available, all have a very small genome size, 
from 25.421 (PorCoV HKU15) to 26.674 (MRCoV HKU18) 
kb, the smallest among all CoVs. Only one papain-like protease 
domain is observed in the nsp3 gene of their genomes. As for 
their gene contents, ORF NS6 was present between the M and N 
genes, and one to four ORFs were also observed downstream to 
the N gene. As for the TRSs, they all have the same putative TRS 
of ACACCA and separation of the TRS from the AUG of the S 
gene by a long stretch of nucleotides. Despite these similar 
genome characteristics among members of Deltacoronavirus, 
NHCoV HKU19 and WiCoV HKU20 possessed genomic fea- 
tures distinct from the other members of Deltacoronavirus, in- 
cluding the amino acids upstream of the putative cleavage sites 
at the junction of nsp2/nsp3, nsp3/nsp4, and nsp4/nsp5. It is 
also notable that NHCoV HKU19, WiCoV HKU20, and 
CMCoV HKU21 occupied the first three branches in the phyloge- 
netic trees constructed using 3CL?”°, Hel, RdRp, and N, indi- 
cating that they could be more ancestral than the other mem- 
bers. Furthermore, these three CoVs were found in large birds, 
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including black-crowned night heron, Eurasian wigeon, and 
common moorhen, in contrast to BuCoV HKU11, ThCoV 
HKU12, MunCoV HKU13, WECoV HKU16, SpCoV HKU17, 
and MRCoV HKU18, which were found in small birds, includ- 
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FIG 4 continued 
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ing bulbuls, blackbird, gray-backed thrush, munias, Japanese 
white-eye, Eurasian tree sparrow, and oriental magpie robin. 
We speculate that the change in genome characteristics (e.g., 
acquisition of s2m) could have occurred during interspecies 
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scale bar. Virus name abbreviations are the same as those in the legends of Fig. 1. 


jumping of the CoV within the large birds before the jump to 
the small birds. Interestingly, the fact that PorCoV HKU15 and 
SpCoV HKU17 are the same species implies that interspecies 
jumping from birds to pigs may have occurred relatively re- 


cently. It is possible that a deletion of 3’ Ns7a and Ns7b had 
occurred during interspecies jumping from birds to pigs, which 
is similar to the observation of interspecies jumping of SARS- 
CoV from civets to humans, with the deletion of 29 bp in ORF 
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8 (25). As for the Asian leopard cat coronavirus, with only the 
Hel, S, E, M, and N gene sequences available, the sequences of 
these gene fragments differ from the corresponding ones in 
PorCoV by less than 2.1% nucleotides or 1.7% amino acids, 
including that for the S gene, which is responsible for receptor 
binding. BEAST analysis showed that the CoV jumped from 
birds to mammals around 523 years ago (Fig. 5). The mixing of 
birds, pigs, and other mammals in domestic environments and 
wildlife markets as well as their close contacts with humans 
may provide the correct environment for interspecies jumping 
and could subsequently pose risks of further genetic changes 
for adapting to human host as in the case of SARS (5). More 
extensive epidemiological studies in different varieties of mam- 
malian species in other parts of the world for members of Del- 
tacoronavirus would further improve our understanding on the 
diversity of this genus as well as its evolutionary history. 
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